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Abstract 

A compelling feature of relativistic mean-field phenomenology has been the 
reproduction of spin-orbit splittings in finite nuclei after fitting only to equi- 
librium properties of infinite nuclear matter. This successful result occurs 
when the velocity dependence of the equivalent central potential that leads 
to saturation arises primarily because of a reduced nucleon effective mass. 
The spin-orbit interaction is then also specified when one works in a four- 
component Dirac framework. Here the nature of the spin-orbit force in more 
general chiral effective field theories of nuclei is examined, with an emphasis 
on the role of the tensor coupling of the isoscalar vector meson [yS) to the 
nucleon. 

PACS number(s): 21.30.-x,12.39.Fe,12.38.Lg,24.85+p 
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I. INTRODUCTION 



The concepts and methods of effective field theory (EFT) have recently elucidated 
the successful nuclear phenomenology of relativistic field theories of hadrons, called quan- 
tum hadrodynamics (QHD) [^|-|^. The EFT framework shows how QHD models can be 
consistent with the symmetries of quantum chromodynamics (QCD) and can be extended 
to accurately reproduce its low-energy features. The EFT perspective accounts for the suc- 
cess of relativistic mean-field models and provides an expansion scheme at the mean-field 
level and for going beyond it []8|,p!0[ . 



One of the most compelling features of QHD mean-field phenomenology is the reproduc- 
tion of spin-orbit splittings in finite nuclei after fitting only to the equilibrium properties 
of infinite nuclear matter. This successful result occurs when the velocity dependence of 
the equivalent central potential that leads to saturation arises primarily because of a re- 
duced nucleon effective mass. The spin-orbit interaction in nuclei is then also specified by 
the structure of the four-component Dirac spinors |p. Does this connection survive in the 
extended EFT framework? Here we re-examine the nature of the spin-orbit force in chiral 
effective field theories of nuclei. 

The strong correlation in QHD models between the effective mass Mq of the nucleon in 
equilibrium nuclear matter and the strength of the spin-orbit force in nuclei has been well 
documented [pAl-[13|,|9[| . Typical values of Mq/M in successful models are roughly 0.6, which 
yield spin-orbit splittings close to experimental values (e.g., 6MeV for p states in ^^O). In 
Ref. Ijl^, a wider class of models, which included variations of the linear sigma model as well 
as all of the standard QHD models used phenomenologically,Q was studied systematically. 
Accurate reproductions of observables in finite nuclei, including the spin-orbit splitting in 
particular, tightly constrained the value of Mq/M to the range 0.58 < Mq/M < 0.64. 

In Ref. [H, an effective hadronic lagrangian consistent with the symmetries of QCD and 
intended for application to finite density systems was constructed. This involved a different 
strategy than in conventional QHD phenomenology. Rather than maximizing predictivity 
by minimizing the free parameters (as in renormalizable models), the goal was to test a 
systematic expansion for low-energy observables, which included the effects of hadron com- 
positeness and the constraints of chiral symmetry. The degrees of freedom are (valence) 
nucleons, pions, and the low-lying non-Goldstone bosons. A scalar-isoscalar field with a 
mass of roughly 500 MeV was also included to simulate the exchange between nucleons of 
two correlated pions in this channel. The lagrangian was expanded in powers of the fields 
and their derivatives, with the terms organized using Georgi's "naive dimensional analysis" 
(NDA) |lT5|,3Jf6|l • The coefficient of each term is written as a combination of the pion-decay 
constant ~ 93 MeV and a larger scale 0.5 ^ A ^ 1 GeV times a dimensionless coupling 
constant. The effective lagrangian is said to be "natural" if these dimensionless coefficients 
are of order unity. 

The result is a faithful representation of low-energy, non-strange QCD, as long as all 



These "conventional" QHD models typically include cubic and quartic self-couplings of the scalar 
field and sometimes a quartic self-coupling of the vector field, and they may also include one-loop 
vacuum corrections. 
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nonredundant terms consistent with symmetries are included. In addition, the mean-field 
framework provides a means of approximately including higher-order many-body and loop 
effects, since the scalar and vector meson fields play the role of auxiliary Kohn-Sham poten- 
tials in relativistic density functional theory . Fits to nuclear properties at the mean- field 
level showed that the effective lagrangian could be truncated at the first few powers of the 
fields and their derivatives, with natural [0(1)] coefficients for each term. Of the two full 
parameter sets identified, one (model G2) provided an excellent fit, including spin-orbit 
splittings, with a nuclear matter value of Mq/M ^ 0.66. More recently, new parameters 
sets with excellent fits have been found with Mq/M ^ 0.70 



An analogous study was made of relativistic "point-coupling" (PC) models in Ref. []T0 
In these models, non-Goldstone mesons are not included explicitly. One instead has an 
expansion of the nucleon scalar and vector potentials in powers and derivatives of local nu- 
cleon scalar and vector densities. As in Ref. |§] , an effective lagrangian consistent with chiral 
symmetry was constructed and fits to nuclear observables were made to test naturalness. 
Excellent fits were found with unprecedented values of Mq/M, as high as 0.74 at nuclear 
matter equilibrium. 

Our purpose here is to examine the nature of the nuclear spin-orbit force in these gener- 
alized models. We find that the most important new feature affecting spin-orbit splittings 



in the lagrangians of Refs. ^ and |]T0[ is the inclusion of a tensor coupling of the omega 
field to the nucleon (or its point-coupling analog). Because the effective lagrangian includes 
a derivative expansion, the tensor term appears at a relatively low order. None of the other 
new terms has a comparable impact on the spin-orbit strength. 

We find a trade-off between the size of the scalar potential (or equivalently, the effective 
mass M*) and the size of the tensor coupling. The existence of such a trade-off has been 
noted in previous studies of the isoscalar tensor coupling in QHD models, but good fits to 
nuclei were found only with relatively small values of the coupling and thus small Mq ^ 
0.6M [^,^. Furthermore, this coupling is usually taken to be zero in one-boson-exchange 
potentials [|1^ or limited to small values due to constraints from free nucleon form factors 
and the assumption of vector meson dominance ||20|-23|. However, as an effective coupling 



in nuclei, which might absorb higher-order many-body effects at the mean-field level, the 
isoscalar tensor coupling could be much larger than these constraints dictate and still be of 
natural size. 

In fits made using more complete meson and point-coupling models, larger (but still 
natural) values are favored, particularly in the point-coupling models |T^. In this work, we 
show that the increase in Mq/M in these models is well accounted for by the contribution 
from the tensor coupling. First, we consider two-component reductions of the models, which 
are reviewed in Sects. and |T|. The spin-orbit potential in the two-component single- 
particle hamiltonian manifests the infiuence of both Mq/M and the tensor coupling, and 
a simple local-density approximation leads to quantitative predictions for the spin-orbit 
strength in nuclei. We find a strong correlation between the size of the isoscalar tensor 
coupling and the increase of Mq/M above 0.6. This discussion applies to a very general class 
of QHD models, including the quark-meson-coupling (QMC) model and the Zimanyi- 
Moszkowski (ZM) model |^ (both of which are found to be deficient). 

A second approach is to transform the nucleon fields in models with tensor couplings. A 
feature of effective lagrangians is the freedom to perform field redefinitions without changing 
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the physics (if all nonredundant terms have been included). It is not possible to completely 
transform away the tensor interaction in a QHD effective lagrangian. To study the con- 
sequences for nucleon-nucleus scattering, however, we can perform an energy- dependent 
transformation to replace the potential in the Dirac equation with an equivalent potential 
(that is, one that generates the same scattering observables) with no tensor term [^. The 
new scalar and vector potentials can then be compared to empirical optical potentials. 

It is well known that in relativistic models applied to hypernuclei, the isoscalar tensor 
coupling of the omega to the hyperon plays a critical role in reducing the strength of the 
hyperon spin-orbit potential p7|P^ . Quark- model arguments imply that the anomalous cou- 
pling is equal to (minus) the vector coupling, which leads to a strong cancellation between 
contributions to the spin-orbit force and therefore small splittings, as needed phenomeno- 
logically. To our knowledge, there is no corresponding implication for the tensor coupling of 
the omega to nucleons.0 



II. SPIN-ORBIT FORCE AND M* WITHOUT TENSOR COUPLING 

In this section, we explore the relationship between the spin-orbit force and the nu- 
cleon potentials in models without tensor potentials. The connection between the spin-orbit 
splittings and the scalar and vector potentials is most easily seen by making a reduction 
to a two-component spin-orbit hamiltonian. We can write the isoscalar part of the Dirac 
single-particle hamiltonian (for spherical nuclei) as 

ho = -iV ■ a + (3(^M -^{r)) +W{r) , (1) 

where M is the nucleon mass, $ is the scalar potential, and W is the vector potential. A 
tensor potential is considered in the next section. Note that we need not assume that $ is 
simply proportional to a scalar meson field 0. In fact, $ could be proportional to (p (as in 
conventional QHD models), or could be expressed as a sum of scalar and vector densities 
(as in relativistic point- couplings models), or could be a nonlinear function of (as in the 
QMC and ZM models). 

To leading-order in 1/M, the spin-orbit term in a Foldy-Wouthuysen reduction of ho 
takes the form 0] 

(t L . (2) 

The expansion is in powers of momenta (i.e., derivatives) and potentials divided by the 
nucleon mass. 

The qualitative impact of the potentials on the spin-orbit force is clear from this expres- 
sion and leads to some basic observations common to all successful models: 



i 

AM^r 



d(P dW 
dr dr 



+ 0(1/M^ 



^We note, however, that in the SU(2) limit of a simple constituent quark model, in which the 
quarks have only Dirac magnetic moments and couple directly to the photon, the nucleon isoscalar 
magnetic moment is small |29|. 
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Successful models are natural p|,[TO[], which implies that the potentials closely follow 



the baryon density. This result follows from the equations for the potentials and the 
close relationship between scalar and vector densities, 

*M^^,.M + o(^.^.^) . (3) 
»'M = yf''»M+o(^.^j , (4) 

where and a„ are 0(1). The precise expressions for the constants and a„ will 
depend on the model; for example, in conventional QHD models as = glf'^/m^ and 
(^v = 5'v/^/'"^v- every case, however, naturalness implies that the corrections are 
order $/M ^ 1/3 smaller than the leading term in the interior of nuclei and smaller 
still on the surface. Furthermore, the scalar density and baryon density are equal to 
within 5 to 10% at ordinary nuclear densities. 

2. Thus a simple local density approximation for $(r) and W{r) should be quite reason- 
able: 

$(r) ^ , (5) 

Po 

W{r) ^ Wo^^ , (6) 
Po 

where $o and Wq are the values of the potentials in nuclear matter at equilibrium 
density po- The remaining factor from the expectation value of /^ls.fw will be universal 
for different models, as long as they accurately reproduce nuclear density profiles. 

3. The Hugenholtz-van Hove theorem relates $o and Wq through the values of the chem- 
ical potential /i and the equilibrium Fermi momentum fcp: 



/i = Wo + V^F + , (7) 

with Mq = M — $0- Since fi and kp are the same in all successful models to very good 
approximation the dependence on Wq can be eliminated in favor of $o- 

4. Therefore, the strength of the spin-orbit interaction and the resulting magnitude 
of the splitting in nuclei should be, to a good approximation, directly proportional to 
the size of $0; the scalar potential in nuclear matter. 

If this analysis is quantitatively valid, we would expect to find that plots of spin-orbit 
splittings in nuclei versus equilibrium Mq/M for different models that reproduce nuclear 
density distributions would lie roughly on a straight line. Such plots are shown in Figs. ^]-^ 
for splittings in ^^O, ^°Ca, and ^ospb using a large number of models. The open circles are for 
QHD-I models PUIH (the original linear Walecka model plus cubic and quartic scalar self- 



interactions). The circles marked with crosses also include one-loop vacuum corrections; this 



is essentially equivalent to adding a particular quintic scalar self- interaction term . Each 
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QHD-I model was fit to a standard set of nuclear matter saturation conditions (equilibrium 
density, binding energy, and symmetry energy) and the charge radius of ^°Ca [^. The 
diamonds are for QHD models from Ref. [^], squares are for point-coupling models from 
Refs. and |T^, and triangles are for QMC and ZM models from Refs. [0] and ||33|. (The 
models with filled symbols contain tensor interactions and are discussed in Sect. JTT[) It is 
evident that there is a universal curve, but it is not a straight line except for Mq very close 
to M (i.e., small in contrast to the conclusion reached in point (4), above. 

For a more quantitative local-density estimate of the spin-orbit splitting, we must go be- 
yond leading order in the Foldy-Wouthuysen reduction, which is not an efficient expansion at 
finite density, even though |p|/M is quite small. The problem is that the Foldy-Wouthuysen 
expansion also involves $/M and W/M, which are less than unity but not negligible. 

Reinhard has proposed a more useful expansion, which counts powers of a (relatively) 
small nucleon velocity v [0. Specifically, 



e„/M = {E^ - M)/M 
($ - W)/M 
M*/M= {M-^)/M 
cr • p M (7 • p 

M M + Ea-<^-W M 



order f ^ , 
order f ^ , 
order 1 , 

order f ^ . 



(8) 



Here Ea is the energy eigenvalue for state a, including the rest mass. Since nuclei are weakly 
bound, the kinetic energy proportional to balances the net potential energy ^ $ — W, 
which leads to the first two results. But while the difference of the potentials is counted as 
a small quantity, the potentials themselves are not. On the other hand, they are still small 
enough compared to M to count M*/M as order 1, which explains the final two counting 
rules. 

Reinhard's procedure constructs order-by-order in v"^ a single-particle hamiltonian /inr 
that acts on two-component (normalized) wavefunctions |Tl|. The latter are related to the 
upper components of the Dirac four-component wavefunctions by a normalization operator. 
He equates the expectation value {ho — E^) calculated with the four-component wavefunc- 
tions with the corresponding expectation value {h^n — ea) calculated with the two-component 
wavefunctions. This condition specifies /inr at each order in w^, up to corrections of the next 
order. 

We have considered two possibilities for the spin-orbit part of h^^ to order f^. 



h 



LS,R 



1 d 



2r dr M{r] 



+ 0{v' 



1 1 dW\ , 

AM r \dr dr I 



cr • L 



(9) 



where M(r) is defined as 



M(r) = M - - ($(r) + Vr(j 



(10) 



and 



hi. 



1 d 1 

2r dr M*(r' 



cr • L 



1 1 rf$ 



2M*2 r dr 



+ 0{v') 



cr • L 



(11) 
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FIG. 1. Spin-orbit splitting for the proton p states in ^^O vs. equilibrium Mg/M for a variety 
of models. The shaded band is an estimate of the experimental uncertainty. Open symbols are for 
models with tensor coupling /v = and filled symbols are for models with /v 7^ 0. The dotted 



curve follows from Eq. (14). See the text for descriptions of the models. 
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FIG. 2. Same as Fig. [l| for the proton d states in ^'^Ca. 
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We have also constructed the 0{v'^) corrections to each. By taking exphcit matrix elements 
of /Ils.r and /i'ls r using two-component wavefunctions for spin-orbit pairs and comparing 
to the energy eigenvalues for the corresponding Dirac wavefunctions, we can determine the 
accuracy of these truncations. We find that /ils,r systematically underpredicts the splittings 
by roughly 5% while /i^s r overpredicts the splittings by about 15%. The predictions using 
the 0{v^) hamiltonians are quite similar in the two cases and are within 3% of the Dirac 
splittings. We will use /iLs.R to 0{v'^) in the sequel. 

We can estimate spin-orbit splittings as follows. We use the local-density expressions 
from Eqs. and (|^) in the numerator of Eq. and replace M in the denominator by an 
average value 

(M) = i(M+(M*)-(iy)) . (12) 

We use an interpolated value of M*: 

(M*) = M* + y{M - M*) , (13) 

and Eq. (0) to define (W) given (M*). One would expect y to be given roughly by ?/ ~ 
1 — {p) / po, where (p) is an average density seen by the nucleons making up the spin-orbit 
pair. In '^''Ca, this would imply y ^ 0.2 for the p states and y ~ 0.4 for the d states. So the 
energy splitting takes the form 



Csplitting = C —=^d{(T • L) , (14) 



where the constant C depends on the average of {l/r)dp^/dr over the nucleus (which is 
essentially the same for the members of a spin-orbit pair), and 5{cr • L) is the difference in 
the expectation values oi cr • L for the spin-orbit pair. The nuclear matter value Mq follows 
from M*/M using Eqs. (0) and (0). 

This local-density estimate of the splitting as a function of nuclear matter Mq/M is 
shown in Figs. as a dotted line. The curves are normalized to the point with the largest 
value of Mq/M to absorb the constant C and any overall normalization errors. The least- 
bound states are well fit with y = 0.4, and the deeply bound states, with y = 0.2. The 
local-density estimates account quantitatively for the spin-orbit trend.0 

The shaded regions in Figs. 0, 0, and ^ indicate the empirical spin-orbit splittings. Since 
we do not calculate rearrangement effects, there are intrinsic uncertainties in comparing 
to experimental splittings.^ We have therefore been conservative in defining these regions. 



^Note that fitting to a different charge radius for ^''Ca shifts the dotted curve up (for a smaller 
radius) or down (for a larger radius). That is why the open diamonds tend to be slightly above 
the curve, since the ^'^Ca radius used is 3.45 fm rather than 3.48 fm as in the other models. 

■^In addition, the Kohn-Sham nature of our calculation implies that the single-particle energy 
eigenvalues do not correspond precisely to observables p. Nevertheless, we expect that differences 
in eigenvalues are good estimates of the splitting, at least near the Fermi surface. 
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However, the intersection of the universal curve with each region imphes that models with 



effective masses outside the range 0.58 < Mq/M < 0.64 proposed in Ref. |14] will not be 
able to reproduce empirical spin-orbit splittings (when there is no tensor coupling). This is 
true of the QMC and ZM |25| models, which are indicated by triangles in Figs. 



One of the appealing features of relativistic models is that one can come to a similar 
conclusion about Mq/M without considering spin-orbit splittings. For example, the study 
by Rufa et al. |]TB[ of conventional QHD-I models (i.e., the linear Walecka model plus cubic 



and quartic scalar meson self-couplings) revealed that a good fit to bulk nuclear properties 
excluding spin-orbit splittings still required Mq/M to be close to 0.6 at equilibrium. 

This result can be understood in terms of a key feature of relativistic models: the ad- 
ditional saturation mechanism that arises from a Lorentz scalar interaction. The velocity 
dependence in the interaction due to the scalar potential appears in the energy of nuclear 
matter (to leading order in k"^) as a modification of the kinetic energy per particle: 



\ \ S k^ f M 



2M \2M* 5 2M\M*(pe 



(15) 



The extra density-dependent repulsion from the M/M* factor leads to saturation 0; in 
successful mean-field models (which include QHD-I models with 0.58 < Mq/M < 0.64) it is 
the dominant saturation mechanism. In these fit to the nuclear matter equilibrium 

density and binding energy alone constrains Mq/M and the implied velocity dependence; 
the four-component structure of the Dirac wave functions then automatically generates an 
appropriate spin-orbit force in finite nuclei. In contrast, in QHD-I models with larger values 
of Mq/M, the empirical equilibrium point is obtained from repulsive three- and four-body 
contributions that become more and more important as Mq/M increases. 

In summary, models without isoscalar tensor couplings exhibit a strong correlation be- 
tween the effective nucleon mass at nuclear matter equilibrium density and the spin-orbit 
splittings in nuclei. Comparison to experiment reinforces the constraint 0.58 < Mq/M < 
0.64 for successful models of this type. Next we consider the generalization that includes 
tensor couplings. 

III. TENSOR COUPLING AND SPIN-ORBIT SPLITTING 

We generalize Hq from Eq. (0) by adding a tensor potential T(r): 

= -iV • CK + /3(m - <l>(r)) + W{r) - i(3oL ■ rT{r) . (16) 

A tensor potential of this form arises in a meson model from the tensor coupling of an 
isoscalar vector meson (a;) to the nucleon. Using the conventions of Ref. [§] , 

Aensor = -J^^a^'N (5^K " d,V^) , (17) 



which implies 



nr) = 4^ . (18) 

^ ' 2M dr ^ ' 
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(Note that /v as defined here corresponds to the combination f^jg 
one-boson-exchange models 
Ref. [|Tg is 



191.) The analog of Eq. ([T 



typically appearing in 
in the point-coupling model of 



/v 



tensor 



which implies 



T(r) 



fv 



(19) 



/2A dr 



(20) 



Thus we can make the rough correspondence /v ^ 2/^. [Note that a spin-2 meson in the 
(1,1) representation is described by a symmetric tensor f^'^. Thus it does not couple to 
Na^^N and will not produce a tensor potential at the mean-field level.] 
We repeat the procedure that led to Eq. (pi) to find 



1 l(d^^dW\^J^l6W^ 



a ■ L 



4M^ r \dr ' dr ) ' 2MM r dr 
and the result for point-coupling models is found from the substitution 

dW 2l dp. 



(21) 



M dr 



A/2 dr 



(22) 



Now suppose that we want to know how M* in a model with fv ^ compares to the value 
Mj^^Q in an "equivalent" model with /v = 0. If we define "equivalent" to mean that they 
have the same spin-orbit strength, we can use Eq. (^1]) to make a rough comparison. Let 
us assume the local-density relations from Eqs. (|^) and (1), so that if we equate /ils.t in the 
two models, the r dependence cancels (in an averaged sense), and we have the condition 



M-{M) 



{MY 



M-{M) , MW) 



/v=0 



{MY 



+ 



M{M) 



(23) 



The analogous condition for the point-coupling models follows after applying Eq. 

We can use Eqs. ( p3D and (^ to transform the region 0.58 < (Mo/M)j^=o ~ 0.64, which 
applies to /v = 0, to the "equivalent" region for > or > 0. Values of (M) and 
(M)/^=o are defined by Eq. (pj) in terms of Mg* and iM^)f,=o using Eq. (|13D with y = 0.2, 
and then (W) and {W)f^=Q follow from Eq. (^. One then inverts the resulting relation for 
Mq given an input (Mo)/^=o- 

The transformed regions are shown in Figs. ^ and along with diamonds and squares 
marking the generalized (/v 7^ 0) models from Refs. 0, ||10| and [|1^. While the shaded 
regions are not definitive indicators of the relation between the tensor coupling and M* in 
successful models, they should provide a reliable estimate of the basic trend. Evidently, the 
bands are consistent with the models with nonzero tensor coupling that were fit directly to 
nuclei. These results are also consistent with new fits in the point-coupling models from 
in which /v is held fixed at zero. The new sets have (Mq/M)j^=o between 0.60 



Ref. \W 



and 0.64. 
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FIG. 5. Correlation between the isoscalar tensor coupling /v and the effective mass in successful 
mean-field models with mesonic degrees of freedom. The shaded region shows the range Mq /M for 
a given value of /v that corresponds to the range 0.58 < Mq/M < 0.64 with = 0, found using 
Eq. (|2^). The QHD models, shown as filled diamonds, are direct fits to nuclear observables from 
Refs. H and 
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FIG. 6. Correlation between the isoscalar tensor coupling and the effective mass in successful 
mean-field point-coupling models. The shaded region shows the range Mg /M for a given value of 
/v that corresponds to the range 0.58 < Mq/M < 0.64 with = 0, found using Eq. (|^). The 
point-coupling models, shown as filled squares, are direct fits to nuclear observables from Ref. pO| . 
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As part of an extensive optimization analysis of conventional QHD models, Rufa et al. 



considered including an isoscalar tensor coupling |T^. They noted the trade-off between 
the magnitude of the spin-orbit splitting and the size of the scalar field. However, their 
chi-squared optimization yielded only a small improvement from the new term,0 and the 
best value had a small value of /v ~ 0.3. (Note that their conventions differ from ours by a 
minus sign.) Reinhard concludes that including the tensor would make a difference only in 
the details, but not in the qualitative picture ITT}]. 



In Ref . , Biro and Zimanyi propose adding a tensor coupling of the u meson to the 



nucleon to the lagrangian proposed by Zimanyi and Moszkowski [^|, in order to bring the 
spin-orbit force up to the strength implied by nuclear spin-orbit splittings. They perform 
a similar expansion to the one above. However, the value of Mq/M in the ZM model is 
roughly 0.85. Figure |] shows that this would require a very large value of /v, which does 
not appear to be compatible with good fits to nuclear observables. 

What constraints exist on the value of /v? In free space there are constraints from nucleon 
electromagnetic form factors under the assumption of vector dominance. Depending on the 



form of vector dominance (i.e., tensor or vector representation |]35|), we obtain a limit on 
/v from the isoscalar anomalous moment of the nucleon or from the isoscalar magnetic 
radius of the nucleon. In either case, one finds the small value |/v| ^ 0.2 [pD|-pB|. From 



Fig. ^ we see that such a limit implies an unimportant change in M* from the results with 
/v = 0. In one-boson-exchange potentials, the isovector tensor coupling is important, but 
the isoscalar tensor coupling is usually taken to be zero [1^. On the other hand, without 



further assumptions, the value of /v in the point-coupling models is not related to the 
electromagnetic structure of the nucleon and is therefore unconstrained. 

It is important to remember, however, that the mean-field parameters obtained from 
fitting to finite nuclei can be very different from the "free-space" tree-level values. These 
mean- field parameters must incorporate the effects of many-body physics beyond the Hartree 
level, which could induce strong renormalizations. We anticipate that isovector couplings 
will be less affected (since the free-space values are already of natural size and nuclei have 
roughly equal numbers of neutrons and protons), and therefore argue that fitting isovector 
parameters using free-space electromagnetic form factors is a good approximation. However, 
the situation with isoscalar parameters is unclear. We might even expect that in the absence 
of symmetry constraints in the medium, they take natural values. This would mean that 
any values shown in Figs. |^ or || are accessible (here, as in Ref. P], /v/4 is the natural 
combination) . 

One might hope to constrain the value of /v from the energy dependence of the optical po- 
tential. A single-particle hamiltonian of the form Hq in Eq. (|Tp implies an energy- dependent 
optical potential for nucleon-nucleus scattering of the form 

Kpt(r; $, W) = 2EW{r) - W{ry + M*{rf - . (24) 

If the potential Kpt('") follows the nuclear density, the slope in energy is given by 2Wo [see 



^The larger values of /v and /v from optimized fits in Refs. P,17,p^ also yielded only a small 



improvement in x (a reduction of about one out of forty) compared to similar models with these 
parameters fixed at zero. 
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Eq. @], and is constrained by low-energy neutron-nucleus scattering to be about 0.6 M 
Il36|,|6[|. To be consistent with experiment, Wq must be as large as the values implied by 
0.58 < Mq/M < 0.64 through Eq. (^. The larger values of Mq, and consequently smaller 
values of Wq, when > (or > 0) could set a limit on the magnitude of the tensor 
coupling. 

We cannot identify directly from hi because of the tensor potential. However, we 
can apply field redefinitions to the nucleon wavef unctions, so that the transformed version 
of Eq. (|l^) does not contain a tensor potential. Such a transformation was applied to Dirac 
wavefunctions in Ref. Given a solution that satisfies hiiplr) = Etplr) for scattering 



states {E > M), we define a transformed wavefunction 0(r) by 

0(r) = e-^(")^»V(r) , (25) 

where U{r) = f^W{r)/2M in the meson model and U{r) = f^p^{r) / f^K in the point- 
coupling model. The function U{r) goes to zero rapidly outside the nucleus, so scattering 
observables are unchanged. 

The new wavefunction satisfies 

-icx-V + We^^^"' + /5(M - $)e2^^°] 0^(r) = Ee^^^Vlr) • (26) 

By rewriting the exponents we can make the identifications: 

M;^^o = cosh 2U + {W - E) sinh 2U , (27) 
Wf^=o = E +{W -E) cosh 2U + M* sinh 2U . (28) 

Therefore the transformation introduces energy dependence into the effective scalar and vec- 
tor potentials. If evaluated in nuclear matter for E ^ fi, the predicted values of {Mq)j and 

(Wo)f=o for the point- coupling models with /v 7^ are consistent with Fig. ^ in particular, 
one obtains 0.58 < {Mq/M)j^^ < 0.61 for each model shown. The predicted values for the 
QHD models in Fig. |^, however, are slightly smaller than expected [(Mq/M)j^=o is roughly 
0.55]. 

If we use Mj^^Q and Wf^=o from Eqs. (pTl) and (pSD to evaluate the optical potential in 
Eq. (g4D, we find 

Kpt(r; $/.=o, Wf^=o) = Kpt(r; $, W) . (29) 



That is, the additional energy dependence in Eqs. (|27D and (|28D cancels out when evaluating 
Vopt- Thus the energy dependence of the optical potential is simply proportional to the 
untrans formed vector potential W and is independent of /v (or /v). If the tensor coupling is 
too large. Figures ^ and ^ imply that Wq will be too small to account for the experimental 
energy dependence. Therefore, smaller values of /v (or f^) are favored, with an upper limit 
around unity. This analysis is not conclusive, however, since the radial dependence of the 
potentials may invalidate the simple argument leading to the comparison with experiment 
HI. This issue will be considered in future work. 
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IV. SUMMARY 



We examined the role of an isoscalar tensor coupling to the nucleon in chiral effective 
lagrangians for nuclei. Effective lagrangians must contain all nonredundant terms consistent 
with the underlying symmetries. Since tensor couplings of the vector mesons to the nucleon 
(or their analogs in point- coupling models) will appear at low order in an expansion in 
powers of fields (or densities) and derivatives, they should be included. 

In models that ignore such couplings, there is a strong correlation between the value of the 
nucleon effective mass Mq at equilibrium nuclear matter density and spin-orbit splittings in 
nuclei. In fact, given a value of Mq, one can predict the calculated splittings of weakly bound 
levels to an accuracy of roughly 0.3 MeV. Moreover, to accurately reproduce the empirical 
splittings, one requires an equilibrium Mq/M between 0.58 and 0.64. Interestingly, similar 
values of Mq are obtained when fits are performed without including information from these 
splittings. 

When an isoscalar tensor coupling is included and allowed to vary freely, fits to nuclear 
properties slightly favor natural-sized values of the coupling, which are significantly larger 
than values deduced from free-space electromagnetic observables. This is particularly true 
in point-coupling models. These larger values of the coupling modify the simple relationship 
between Mq and the spin-orbit splittings noted above; however, it is possible to estimate 
the contribution of the tensor interaction to the spin-orbit force and to show that the larger 
values of Mq obtained in nuclear matter are consistent with the observed splittings. 

We emphasize that the small size of the isoscalar tensor coupling deduced from free-space 
electromagnetic data is not incompatible with the larger values obtained in fits to nuclei, 
since the in-medium parameters must contain many-body effects. In all of the models dis- 
cussed here, the resulting tensor coupling is still natural. It remains to be seen whether 
explicit calculation of the many-body contributions can account for the larger values fa- 
vored by nuclei. Moreover, it may be possible to constrain the tensor coupling by including 
additional observables in the fits; the energy dependence of the optical potential for nucleon- 
nucleus scattering is a promising candidate. 
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